Observation of anomalous spin-state segregation in a trapped ultra-cold vapor 
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We observe counter-intuitive spin segregation in an inhomogeneous sample of ultra-cold, non- 
condensed Rubidium atoms in a magnetic trap. We use spatially selective microwave spectroscopy 
to verify a model that accounts for the differential forces on two internal spin states. In any simple 
understanding of the cloud dynamics, the forces are far too small to account for the dramatic 
transient spin polarizations observed. The underlying mechanism remains to be elucidated. 
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Imagine a near-ideal gas consisting of a mixture of 
two distinct molecular species. Apply a differential force 
that tends to pull the two species apart. What de- 
termines whether the two species segregate, or remain 
mixed? From a kinetics point of view, the question be- 
comes whether the segregating drift velocity, a function 
of the differential force and the mean free path, is larger 
or smaller than the effective velocity of remixing diffu- 
sion. The "diffusive velocity" depends on the character- 
istic size of the sample. The equilibrium point of view 
focusses instead on whether the energetic benefit derived 
from separating in the differential potential outweighs the 
corresponding entropic cost. Again, for a given magni- 
tude force, the ultimate outcome hinges on the size of 
the sample. Both points of view lead to the same quan- 
titative criterion - is the differential potential across the 
sample large or small compared to the thermal energy 
ksT? As an example, a mixture of gaseous oxygen and 
nitrogen in a room-sized box does not undergo separation 
under the differential force of Earth's gravity, because the 
gravitational potential of the differential mass (4 amu) , 
evaluated over a few meters, is very small compared to 
k B T. 

In this paper, we study the segregation behavior of 
a gaseous mixture. Our system consists of magneti- 
cally trapped ultra-cold rubidium gas in two distinct spin 
states. These states can be interconverted in the pres- 
ence of an applied microwave magnetic field, but in the 
absence of the applied drive the two states do not in- 
tercom/erf and can be thought of as two distinct fluids 
K|. Inter-atomic interactions and differential magnetic 
moments give rise to a small differential potential that 
would tend to separate the two species, except that it is 
a thousand times smaller than ksT. In contrast to the 
above arguments we observe that the initially homoge- 
nized sample nonetheless undergoes a transient segrega- 
tion. 

In the first part of this paper, we describe the use of 
spatially selective Ramsey spectroscopy to characterize a 
differential potential acting on the two spin states. In the 
second part, we describe our observations of the atomic 



motion induced by the differential potential over longer 
timescales. The atomic motions are much too large to 
be described simply as arising from an acceleration from 
the differential potential. 

The experimental apparatus includes a single magneto- 
optical vapor-cell trap which precools a sample of 87 Rb 
atoms for transfer into a quadrupole magnetic trap. The 
electro-magnetic quadrupole trap coils are mounted on a 
servo-controlled linear track that mechanically moves the 
atoms to a UHV region where they are transferred into 
a Ioffe-Pritchard magnetic trap for further cooling via 
rf evaporation. Our Ioffe-Pritchard trap uses permanent 
magnets to produce strong radial confinement. Electro- 
magnetic coils generate a stable and adjustable bias field 
and a weak confining potential in the axial direction. The 
trapping frequencies are {7, 230, 230} Hz. Typical sample 
sizes are 10 6 atoms at a temperature of 600 nK (factor of 
2 above the Bose-Einstein condensation temperature). 
The two hyperfine states of interest are |1) = \F = 
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-1) and |2) = \F = 2, 
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We coherently 



couple the two hyperfine states via a two-photon drive 
P). The drive consists of a microwave photon (~ 6.8 
GHz) and a rf photon (1-3 MHz). The sum frequency is 
detuned 200 Hz from the transition, and the effective two- 
photon Rabi frequency is 2.5 kHz. The drive is turned 
on for ~ 100 [is to transfer half of the atoms (initially all 
in the |1) state) to the |2) state (a ^ pulse). A second | 
pulse is applied after a variable delay . The final popu- 
lation and spatial distribution in either spin state can be 
separately determined using absorption imaging |J. 

We determine the transition frequency between states 
|1) and |2) by varying, in each realization of the experi- 
ment, the delay time between the two ^ pulses and mea- 
suring the relative population in the |1) state. Typical 
delay times range from 1-60 ms. The frequency of the 
resulting Ramsey fringes is the difference between the 
atomic transition and the oscillator frequency of our two- 
photon drive pj. In order to study the spatial variation 
of the transition frequency between the two states, we 
break up the images into 9 bins along the axial direction 
and integrate the number of atoms in each bin (Fig. [j]). 
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There are two different mechanisms that shift the tran- 
sition frequency. The first, dominant in the limit of low 
atomic density, is the differential Zeeman effect. The 
atoms are confined in the parabolic minimum of the mag- 
netic field and each spin state individually experiences 
a Zeeman shift, which varies by 30 kHz from the cen- 
ter of the cloud to the edge of the cloud ||. States |1) 
and 1 2} were selected because the Zeeman shifts cancel to 
great accuracy. The Breit-Rabi equation H predicts that 
the frequency splitting is first-order independent of field 
at 3.23 G, increasing quadratically with magnetic field 
about that minimum. The magnetic field in turn has a 
minimum value, B^ias, at the center of the magnetic trap, 
and increases quadratically with displacement along the 
trap axis. Thus for a Bbi as = 3.23 G, the differential Zee- 
man shift across the cloud is quartic in axial position; the 
quartic shift from trap center to sample edge is 0.7 Hz. 
For \Btias — 3.23G| > 0.02 G, the dominant field depen- 
dence with position is quadratic. For instance, for B^ias 
= 3G, the frequency shift from cloud center to edge is -7 
Hz. 

The second frequency shift is due to the mean-field 
energy that arises from the self- interaction of the atoms. 
The two states have slightly different s-wave scattering 
lengths and therefore different interaction energies. The 
shift of the transition energy is given by 

AE m t = [2n 2 a 2 2 ~ 2nia n + 2(m - n 2 )ai 2 ] (1) 

m 

where a.^ is the s-wave scattering length between states 
\i) and rii is the density in state \i) and the pref- 
actors of 2 account for the enhancement of s-wave col- 
lisions for non-condensed bosons. For 87 Rb the s-wave 
scattering lengths are an = 100. 9ao, 022 = 95.6ao, and 
ai2 = 98.2ao 0| where ao is the Bohr radius. The mean- 
field shift has very little dependence on relative popu- 
lations in the two states because a\ 2 ~ (an + a 22 )/2. 
After a ^ pulse, ri\ w n 2 w n/2, and to sufficient accu- 
racy AE rn f = (a 22 — aii)A-KTin/m. 

We measure the mean-field shift spectroscopically. The 
sample density is inhomogeneous, with a gaussian profile 
in each dimension. The radial oscillation period is short 
compared to the duration of the Ramsey measurement 
and to the mean collision time (~ 10 ms). Thus the 
random thermal motion of the atoms effectively averages 
over the radial density inhomogeneity. The axial oscil- 
lation period, on the other hand, is long compared to 
the measurement duration and to the collision time. To 
a reasonable approximation one can treat the atoms as 
"frozen" in the axial direction. To further reduce the ef- 
fect of axial density inhomogeneity, we use the data only 
from the center bin, a region of relative uniformity. We 
create clouds with various peak densities and measure the 
transition frequency of the center bin. The slope of the 
frequency vs. density fit (—0.52 x 10 12 Hz/cm 3 ) is more 
negative than a prediction based on photo-associative 



data by 20% 0. However, it is entirely possible that 
our analysis of our cloud images could have systemat- 
ically underestimated the cloud density by as much as 
20%. More troubling is that the zero-density extrapo- 
lation, at 6834678118.0(7)Hz at 3.23 G, is too high by 
4.4(7) Hz. We are unable to account for this systematic 
offset. 

For a particular peak density we measure the axial spa- 
tial variation in frequency by binning the cloud into sev- 
eral sections and determining a frequency for each bin. 
This measurement is done at B^ias = 3.23 G to minimize 
the magnetic contribution to the inhomogeneity. The 
transition frequency vs. axial position along the cloud is 
shown in Fig. 0c. 
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FIG. 1. (a) Images of F=l atoms (upper) and F—2 atoms 
(lower) after two |- pulses. The delay between the pulses is 
listed at the top. The edges of the cloud have become out of 
phase with the center. The black lines show how the cloud 
is divided up along the axial direction to measure the spa- 
tial dependence of the transition frequency. Each image is 1 
mm on the horizontal (axial) dimension, (b) The center-bin 
|1) optical density as a function of time between the two \ 
pulses, (c) The plot shows the spatial dependence of the tran- 
sition frequency in the axial direction. Offset frequency is 
defined as the transition frequency - 6834682610.9 Hz. The 
solid line is a gaussian fit to the data. These data were taken 
at Btias = 3.23 G. 

We can change the frequency shift across the cloud 
and nearly cancel both of the frequency shifts by oper- 
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ating at different a value of Bu as . The cancellation of 
the mean-field and magnetic frequency shifts can not be 
perfect because the cloud has a gaussian density profile 
and the magnetic field shift is essentially parabolic across 
the cloud. However, for a cloud with 1 x 10 6 atoms at a 
temperature of 500 nK there will be less than 1 Hz shift 
across the cloud at B bias = 2.83 G (Fig. ||). 
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FIG. 2. The transition frequencies at different positions 
along the cloud are plotted as a function of Bbias- The cloud 
has a peak density of 4.0 xl0 13 cm~ 3 . Offset frequency is 
defined as the transition frequency - 6834682610.9 Hz. The 
frequency shift at 3.23 G (negligible magnetic contribution to 
the inhomogeneity) has a gaussian shape. At B^ias = 3.37 G 
the magnetic shift and the mean-field shift are in the same di- 
rection and therefore the inhomogeneity is across the cloud is 
increased. When the cloud is in a field of 2.83 G inhomogene- 
ity due to the mean-field shift is almost completely cancelled 
by the magnetic shift. Data points are connected to guide the 
eye. 

Conversely, if we operate at Bu as > 3.23 G, the curva- 
ture of the magnetic contribution to the frequency inho- 
mogeneity adds with the same sign to that from the den- 
sity contribution, and we see enhanced total frequency 
shift across the cloud (Fig. ||). Near the center of the 
cloud, the frequency shift is parabolic with displacement 
from the center, and the curvature of the differential po- 
tential can be conveniently characterized by 



vdijf = l/2-K-y/h/m x (d 2 v 12 /dz 2 ) (2) 

where m is the mass of 87 Rb and 1/12 is the transition fre- 
quency H . Vdif f may be thought of as the axial motional 
frequency of a hypothetical test atom, in state |2), under 
the influence of only the differential magnetic and mean- 
field forces. For the data shown in Fig. 0, measured j j 
ranges from 0.15 Hz for Bi, ias — 3.37 G, down through 
O.li Hz for B Uas = 2.83 G. 

We have used spectroscopic methods to characterize 
the spatial dependence of the differential potential expe- 
rienced by states |1) and |2) in our trap. The data are 



quantitatively accounted for by a simple model involving 
magnetic and mean-field potentials. We now describe the 
observed mechanical effects that the gradients of these 
potentials have on the distribution of the relative atomic 
densities of states |1) and |2). 

In the twin-pulse spectroscopy described so far, the 
first pulse sets up a coherent superposition of the two 
states; the delay time between the two pulses allows the 
relative internal phase, <j), in the coherent superposition 
(|1) + e^|2)) to evolve under the inhomogeneous poten- 
tial, and the second pulse then converts the spatial pat- 
tern in <p into a pattern of relative density for ready imag- 
ing. In the experiments described below, the second J 
pulse is omitted, and the delay between the first pulse 
and the density imaging is increased to allow time for 
the atoms to move under the influence of the develop- 
ing spatial gradient in <f>. The patterns in the density 
distribution which we image must result from the physi- 
cal displacement of atoms during the period after the ^ 
pulse ||. 

The results of this class of experiment are summarized 
in Fig. |j. Starting around 50 ms after the ? pulse, 
the atoms in the |1) state are observed to move outward 
along the trap axis, while the atoms in the |2) state move 
inward. Maximum separation is realized at about 100 
ms, then the two clouds relax back to their original, well- 
overlapped distribution by around 200 ms. No further 
evolution is observed. The spin segregation is dependent 
on the curvature of the differential potential. For Vdiff < 
0.05 Hz (Fig. D,(ii)), no separation is observed above 
the imaging noise. For Vdiff > 0.18 Hz (Fig. ||,(iv)), 
we observe what we call "higher-order" effects, with the 
relative spin-density developing multiple stripes. 
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FIG. 3. (a) Spin segregation appears in the spatial dis- 
tribution of atoms in the |1) state at various times after the 
■| pulse. |1) state atoms move outward, then relax back to 
their initial condition. Initial cloud conditions are peak to- 
tal density 1.8 x 10 13 cm -3 , T=630 nK, v diff = 0.1 Hz. (b) 
Same as for (a), but showing the |2) state distribution. |1) 
and |2)-state motions are complementary, and leave the total 
density profile unchanged, (c) |1)- state axial density profile 
imaged at times after the ^ pulse, listed on the left column 
. For all columns, T = 600 nK. Vdiff values are listed at 
the bottom of each column. The first column of cross sec- 
tions has an imaginary Vdiff corresponding to negative cur- 
vature of the frequency. Initial density (in 10 13 cm -3 ) and 
Bbias (in Gauss) are for each column (i) 1.82, 2.99 (ii) 1.82, 
3.12 (iii)1.82, 3.52 (iv)1.82, 3.79 (v)0.254, 3.35 (vi)1.16, 3.31 
(vii)2.05, 3.27 Columns (i)-(iv) are at fixed density and pro- 
gressively higher Vdiff, showing turn-on of effect with increas- 
ing Vdiff - Columns (v)-(vii) are at fixed Vdiff and at progres- 
sively higher density, showing turn-on of segregation effect 
with density. 

We find the segregation to be startling. While the 
mean-field and magnetic forces can combine to provide a 
small buoyant force that tends to push the |1) atoms out 
towards the lower density regions of the clouds, the mag- 
nitude of the observed effect is entirely disproportionate. 
For Vdiff = 0.1 Hz, for instance, the total frequency shift 
A^i2 from center to edge of cloud is 6 Hz. This should 
be compared to ksT/h ~ 10 kHz. The criterion for sep- 
aration discussed in the introduction is thus clearly not 
satisfied. This objection would be surmounted if the ini- 
tially well-homogenized and non-condensed populations 
in each spin-state nonetheless moved collectively (the rel- 
evant inequality would then be N a t oms (Ai>i2h) 3> kT 
rather than Av^h <C kT) |10|. One would still need 
to account, however, for the rate of segregation. If the 
atoms move ballistically in the differential potential, the 
minimum time for the peak segregation to develop should 
be at least 1/(4^//) [|l(|. But the observed segregation 
times are at least 10 times faster. 

As seen in Fig. ||c, to induce segregation one needs not 
only a sufficiently large value of fdi f / but independently 
one needs adequate density. It is possibly significant that 
the segregation seems to disappear when the density gets 
so low that the cloud is no longer collisionally thick in 
the axial direction. While the observed segregation rate 
is fast compared to 4 Vdiff, it is slow compared to Av z , 
where v z is the axial trap frequency. That is, the random 
thermal speed of the atoms, if it were selectively directed, 
is more than adequate to move the atoms the required 
distance. We speculate that as the internal phase (j> de- 
velops a spatial gradient, it somehow acts as a "switch" 
to control velocity-changing collisions such that atoms in 
the |2) state are preferentially scattered inward and |1) 
atoms outward. 

Finally, the observation of non-diffusive evolution of 
spin populations compels a comparison with the spin- 



wave phenomenon observed |TlJ and discussed exten- 
sively |l^] in the 1980s. We understand the earlier work 
to be an effect of the transverse spin. The present work 
describes evolution of the relative population of energy 
eigenstates, a longitudinal spin effect. Additional study 
may lead to connections between the two spin effects. 

We acknowledge useful conversations with the other 
members of the JILA BEC collaboration and with Kurt 
Gibble. This work is supported by the NSF and by NIST. 
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